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INSTABILITY OF THE 13.8 DAY PERIOD IN X-RAY EMISSION FROM 4U 1700-377 
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ABSTRACT 

We present a new result on long-term periodicity searches of the X-ray emission from 4U 1700-377 
using RXTE ASM and CGRO BATSE data. The hard X-ray data (20 - 200 keV) from early BATSE 
observations (1780 days from JD 2448370, before RXTE observations started) show evidence of a 13.8 day 
periodicity. The long term periodicity became substantially less prominent in the data from later BATSE 
observations and in the soft X-ray data from the ASM (2 - 10 keV) observations. We demonstrate that 
disk precession models can explain the 13.8 day period and its instability in the X-ray emission from 
4U 1700-377. 

Subject headings: X-ray binary, long term period, disk accretion 



1. INTRODUCTION 

First discovered by Jones et al. (1973), 4U 1700-377 is a 
3.4 day X-ray eclipsing binary system. Its optical counter- 
part, HD 153919, is a bright early-type star. The system 
has been extensively studied at both optical and X-ray 
wavelengths. While it is considered to be a neutron star 
binary system there has been no confirmed detection of 
pulsation. 

Konig & Maisack (1997) reported a possible 13.8 day pe- 
riodicity in the X-ray emission from 4U 1700-377. Their 
results are based on the long term observations by the 
All-Sky Monitor (ASM, 2-10 keV) (Levin et al. 1996) on 
the Rossi X-ray Timing Explorer {RXTE) and the Burst 
And Transient Source Experiment (BATSE, 20 - 200 keV) 
(Fishman et al. 1989; Harmon et al. 2002) on the Compton 
Gamma Ray Observatory (CGRO). More data has been 
collected since their report. Here we present the results of 
long-term period searches using all the data available as 
of this writing. 

2. DATA SELECTION AND ANALYSIS 

We employed All-Sky monitor X-ray light curves from 
BATSE on CGRO satellite and ASM on the RXTE satel- 
lite. All data presented here are publicly available by an 
on-line service from the Laboratory for High Energy As- 
trophysics (LHEA) at NASA/GSFC 3 . 

Using Earth occultations BATSE provides long term 
flux histories for bright X-ray sources in the energy range 
from 20 to 300 keV. In this paper, we used 3321 days 
(from JD 2448370) of daily-averaged, long-term flux from 
4U 1700-377 given by the BATSE occultation technique 
(Figure 1). The ASM on the RXTE satellite usually dwells 
on a bright source a few times a day and it covers an en- 
ergy range from 2 to 10 keV (Figure 1). We used 2221 
days (from JD 2450088) of ASM X-ray time series data. 
No additional cut was applied to the data for analysis. 



For the period search, Konig & Maisack (1997) used the 
epoch-folding technique, which is based on Leahy et al. 
(1983). The epoch-folding technique in Leahy et al. (1983) 
does not take into account phase selection effects. Here 
we used the phase-averaged epoch-folding technique (Col- 
lura et al. 1987). The phase-averaged epoch-folding tech- 
nique is less vulnerable to false alarms due to random fluc- 
tuations than a single-phase epoch-folding technique. In 
practice, the outcome of the phase-averaged epoch-folding 
technique is similar to that of the epoch-folding technique 
with an optimally chosen bin size (Collura et al. 1987). 

Figure 2 shows the phase-averaged epoch folding results 
for the light curves of 4U 1700-377 in Figure 1 (15 phase 
bins). In Table 1, we estimated the orbital period (Pq), as- 
suming the sinusoidal signal of light curves (Leahy 1987). 
It should be noted that the true orbital period varies as 
P/P = -3.3 x 10~ 6 /yr (Rubin et al. 1996). 

In comparison with the previous report (Konig & Mai- 
sack 1997), the phase-averaged epoch folding on the 1780 
days of BATSE data (Figure 3) produces roughly the same 
level of the noise fluctuations as in Figure lb of Konig & 
Maisack (1997), but Xred & t the multiples of the orbital pe- 
riod is about twice as large as those of Konig & Maisack 
(1997). This is due to the superior performance of phase- 
averaging in the epoch-folding technique. In the case of 
the first 573 days of the ASM data, the noise fluctuation 
of Xred by the phase-averaged epoch folding is dramati- 
cally smaller than that of the previous reports (Figure 3 
in Konig & Maisack (1997)). In Figure 3, the second har- 
monics of the 13.8 day period is also prominent as a broad 
peak around 27.6 days. 

The first two plots in Figure 4 show a detailed view of 
the P 2 and P 3 regions (P l = (i + 1)P ) from the BATSE 
data. Figure 4 also shows the P3 region of y^ed fr° m the 
first 1780 days and the next 1541 days of the BATSE data. 
The data were divided in this way in order to make a direct 
comparison with the results in Konig & Maisack (1997). 
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Fig. 1.— X-ray light curve from 4U 1700-377 : CGRO BATSE light curve by occultation technique (top : 20 - 200 keV) and RXTE ASM 
light curve (bottom : 2 - 10 keV). The intervals covered in Konig & Maisack (1997) are marked by a horizontal bar. 



Table 1 

Estimate of the orbital period P and its width AP in the epoch-folding diagram 



Data set 


Po estimate 3. 


APo cstimatc a /theory b 


Xrcd 


days 


day 


xKT 3 day 




BATSE 3321 


3.41155(10) 


3.11(5)/3.50 


58.2 


ASM 2221 


3.41161(05) 


5.41(l)/5.24 


366.4 



a P and AP is calculated for the sinusoidal signal (Leahy 1987). 
b AP = P 2 /T, where T is the time span of the data set. 
b The maximum value of the peak 
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Fig. 2.— The phase- averaged epoch folding results of the light curve from the CGRO BATSE data (top: 3321 days from JD 2448370.0) 
and the RXTE ASM data (bottom: 2221 days from JD 2450088.6) (15 phase bins) 
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Fig. 3. — The phase-averaged epoch folding results for the first 1780 days (from JD 2448370) of the BATSE data: compare this plot with 
the figure lb in Konig & Maisack (1997), where X^ eA at Pj is not always larger than that at Pj+i due to the noise fluctuations. One can also 
notice a broad peak around 27.6 days, which is the second harmonics of the 13.8 day period. 
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In Figure 4, one can notice that there is a bump around 
13.8 days in the P3 region (b), which indicates the pres- 
ence of 13.8 day periodicity. Judging from the P2 region 
(a) and the expected width of P 3 , the bump is not from 
the side- lobe of P 3 . The bump is more significant from the 
first 1780 days of data (c) and its peak is at 13.82(2) days, 
which is consistent with the result in Konig & Maisack 
(1997). But there is no clear indication of the 13.8 day 
periodicity in the rest of the data (d) , although the epoch- 
folding search shows some features around 13.8 days. 

Figure 5 shows the P 2 and P 3 regions of the epoch- 
folding results using the 2221 days and the first 573 days 
of the data from the RXTE ASM observations. 

As shown in Figure 5 (b), there is no clear indication 
of the 13.8 day peridocity from the RXTE data. The 
peak around 13.8 days appears to be a side-lobe of P3, 
since the similar feature is present in the P2 regime (a). 
The two bottom plots employ the same 573 days of the 
data as in Konig & Maisack (1997). Judging from the 
width of the P3 peak in the plot and its expectation 
(AP 3 = 4P 2 /T tot w 0.08 day where T tot = 573 days) we 
can conclude that the bump around 13.8 days is largely 
due to the side-lobe of P3. One could argue that the right 
side-lobe near 13.8 days is larger than the left side-lobe 
at 13.5 days in the bottom-right plot, which may be in- 
terpreted as the presence of an additional bump, but the 
similar trend is noticed in other Pj regions. 

We divided the data into several sets (for example, di- 
viding into 700 days of data for each set provides 5 sets of 
BATSE and 3 sets of RXTE data) and applied the epoch- 
folding technique (in the case of RXTE, a separate light 
curve is available for three energy bands). All the results 
(not shown here) indicate that the 13.8 day periodicity in 
the BATSE data is significant only in the first <~ 1700 days 
of the data, which were taken before the RXTE observa- 
tion on the source started. 

Figure 6 shows the folded light curve at the 13.8 day 
period, using the first 1780 days of the CGRO BATSE 
and the first 573 days of the RXTE ASM data. In the 
case of the ASM data, there is a series of small peaks 
caused by the proximity of P3 to the 13.8 day period for 
the total time span (573 days). However, the peak at the 
0.8 phase in the ASM data is similiar to the one in the 
BATSE data, which is considered the main source of the 
high Xred va hie at the 13.8 day period in the epoch folding 
analysis. This implies that the ASM data may contain a 
long term periodicity and its signature is too weak to be 
directly detected, particularly due to the proximity of P3. 

In order to detect the 13.8 day periodicity unambigu- 
ously, it is important to remove the contribution from 
the side lobe of P3 in calculating x^cd f° r a gi ven period. 
Therefore, we applied the epoch folding procedure only 
on the uneclipsed data, the orbital phase of which falls in 
between 0.2 and 0.8. If a folded phase bin at a certain 
period (usually around the multiple of the orbital period) 
contains no data due to exclusion of the eclipsed data, we 
simply ignore the bin for calculation of the x? e d at tnc 
period 4 . 



The result is shown in Figure 7 and the estimations of 
the long term periodicity are summarized in Table 2. The 
figure indicates that the 13.8 day periodicity became weak 
but it is clearly present in the RXTE data and the later 
part of the BATSE data. The peak position of the \red 
around 13.8 day down-shifted and the peak appears to 
have spread, indicating that the peak consists of more than 
one periodicity (Tabic 2). In Figure 7, one can also notice 
its second harmonics as well as the 13.8 day periodicity. 
In the next section, we will discuss the origin of the long 
term periodicity and its weakening. 

3. DISCUSSION 

Konig & Maisack (1997) argued that the origin of the 
13.8 day periodicity is a free precession of the neutron star 
arising from misalignment of the orbital axis and the neu- 
tron star spin axis. The analysis of the CGRO BATSE 
and RXTE ASM data indicates that the 13.8 day period- 
icity weakened dramatically after ~ JD 2450000. If the 
precession of the neutron star caused the 13.8 day peri- 
odicity, the weakening of the periodicity implies that the 
neutron star precession has been dissipating. Such an in- 
terpretation may not be the most favorable one because 
it suggests that 4U 1700-377 was experiencing a sudden, 
unusual change of the system geometry. 

In this section, we investigate accretion disk precession 
models in order to explain the 13.8 day periodicity and the 
instability of X-ray emission from 4U 1700-377, without 
invoking a dramatic change of the system geometry. So far 
about 15 X-ray binary systems are known to show a long 
term period, often called a superorbital period, in their 
light curve, and in many cases, the superorbital period is 
successfully explained by a precession of the accretion disk 
with respect to the orbital plane (Wijers & Pringle 1999; 
Ogilvie & Dubus 2001). Although 4U 1700-377 is a pri- 
marily stellar wind driven system, there are several obser- 
vational results supporting a possible transient accretion 
disk in the system. Therefore we first explore the clues 
and the characteristics of an accretion disk in the system. 
Second, we apply two different accretion-disk precession 
models - radiation driven (Ogilvie & Dubus 2001 - here- 
after OD), and magnetically driven (Lai 1999) warped disk 
precession models and we study the implications. 

3.1. A long-term, transient accretion disk? 

Disk accretion is one of the most radiation-efficient ac- 
cretion processes. When a stellar wind driven system 
like 4U 1700-377 forms an accretion disk, the disk accre- 
tion usually enhances the overall accretion rate, leading to 
higher X-ray emission in general. So a sudden increase of 
X-ray flux from some systems indicates the formation of 
an accretion disk in the system. For example, GX 301-2 
is a wind-driven high-mass X-ray binary (HMXRB) with 
a partially underfilled Roche lobe, similar to 4U 1700-377. 
The observed spin-up episodes in GX 301-2 are success- 
fully explained by the formation of a transient disk and its 
angular momentum exchange with the neutron star (Koh 



In the epoch folding procedure, \ 2 = E(xj — x) 2 /a?, where X{ is the counting rate in the ith phase bin, x the average counting rate. Under 
the null hypothesis, the standard deviations Uj should be defined as erf = <rj? t /raj, where utot is the standard deviation of the unfolded time 
series and rij is the number of data points in the ith phase bin (Larsson 1996). Therefore, when n, = 0, the ith phase bin contributes nothing 
to x 2 . 
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Fig. 4. — The detailed view of epoch folding X^ed from the BATSE data: P2 (a) and P3 (b) region from all the data (3321 days from JD 
2448370.0), P3 region from the first 1780 days (c) and from the next 1541 days (d) of the data. 



Table 2 

Estimate of the long term period P) and its width AP) in the epoch-folding diagram 



Data set 
days 


Pi estimate 
day 


APi estimate/theory 
xlO" 1 day 




BATSE first 1780 


13.82(2) 


0.97(11)/1.07 


4.93 


a first 1780 


13.82(2) 


0.83(18)/1.07 


5.05 


a ' b ncxt 1541 


13.71(3) 


2.24(34)/1.22 


2.80 


3321 


13.81(1) 


1.21(36)/0.57 


3.98 


a < b 3321 


13.81(1) 


1.34(33)/0.57 


4.14 


ASM a 2221 


13.74(1) 


2.60(13)/0.85 


13.5 



a the uncclipsed data (the orbital phase: 0.2 - 0.8) 

b exhibits multiple peaks around 13.8 days. Here we consider them as a part of a single broad peak to estimate the overall 
peak location and the width. 
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Fig. 5. — The detailed view of epoch folding X^ e( j from the ASM data: P2 (a) and P3 (b) region from all the data (2221 days from JD 
2450088.6), and the same from the first 573 days (c,d) of the data 
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Fig. 6.— Folded pulse profiles at 13.8 day; the first 1780 days of the CGRO BATSE data (top), and the first 573 days of the RXTE ASM 
data (bottom), zero is JD 2448369.0 




Fig. 7. — The phase- averaged epoch-folding results on the uneclipsed data: the first 1780 days of the BATSE data (top), the rest of the 
BATSE data (middle), all the ASM data (bottom). The right panels are a close-up around 13.8 days for each case. The long term periodicity 
and its second harmonics are clearly noticeable. 
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et al. 1997). Since the spin-up episodes in GX 301-2 are 
accompanied by enhancement in the X-ray flux from the 
system, an increase in X-ray emission can be used as indi- 
rect evidence of an accretion disk in the system. 

In the case of 4U 1700-377, the intensity distribution of 
the first 1780 days of the BATSE data shows a clear differ- 
ence from that of the next 1541 days of the data (Fig. 8). 
In the first 1780 days of the data, which exhibit the 13.8 
day periodicity rather strongly, there are relatively more 
occurrences of high flux events than in the next 1541 days 
of the data. The difference in the average flux of these 
two data sets is 0.0083 in the given unit. Randomly di- 
viding the combined data into two sets of data - 1780 and 
1541 day long each, we found that the probability of get- 
ting such a difference from random fluctuations in a single 
distribution is negligible (< 10~ 8 ). 

The increase in the X-ray flux may be caused by changes 
in the stellar wind of the companion. However, the corre- 
lation between the frequency of high X-ray flux episodes 
and the strength of the 13.8 day periodicity indicates that 
both may be related to disk accretion. If the difference in 
the average flux is caused by changes in disk accretion, the 
estimated fraction of the X-ray flux from disk accretion is 
> 14%. This is a very crude estimation due to the limited 
data statistics. 

Other compelling evidence of disk accretion is quasi- 
periodic oscillations (QPOs) of X-ray flux. It is often 
believed that QPOs are caused by the material orbiting 
near the inner edge of the accretion disk. In 4U 1700-377, 
QPOs were observed on several occasions during the un- 
cclipsed phase and they lasted ~ 200 - 1200 seconds (Ru- 
bin et al. (1996) and references therein). In order to ex- 
plain these QPO episodes, Rubin et al. (1996) proposed 
a beat frequency model (BFM) mechanism operating in 
a transient accretion disk near the X-ray source. In the 
BFM, a scries of clumps in the accretion disk orbit the 
neutron star at the magnetopause in Keplerian motion, 
and the QPOs are a modulation at the beat frequency 
between the Keplerian rotation frequency at the magneto- 
sphere and the neutron star spin frequency. 

The episodic detection of QPOs may imply a transient 
accretion disk (Rubin et al. 1996). Rubin et al. (1996) es- 
timated 0.8 < il < 1.0 (Roche- lobe filling factor), which 
also supports the idea of accretion through a transient 
disk. Taking the idea of a transient disk a step further, 
we suggest that a transient disk in 4U 1700-377 can be 
quite persistent (> 1000 days). 

First, it should be noted that the episodic behaviour 
could be simply due to the occasional formation of cer- 
tain clumps in a persistent or transient disk. Second, if 
the companion O star is a just-underfilled- Roche lobe sys- 
tem (0 < 1 - fi C 1.0), the star may frequently undergo 
at least partial Roche-lobe overfilling due to various in- 
stabilities. For instance, tidal force driven disruption of 
the atmosphere of the companion star can cause instabil- 
ity. Naturally the time scale for such overfilling would be 
related to the orbital period - 3.4 days. Therefore, such 
instabilities provide a semi-periodic source to feed a tran- 
sient or permanent accretion disk. 

The continual Roche-lobe overfilling may result in ac- 
cretion through several successive rings rather than a con- 
tinuous disk. However, as often assumed for many disk 



models, it is believed that the inward drift rate of such a 
ring is much slower than the rate of viscous-driven ring 
broadening. The material is likely to be accreted by one 
continuous disk rather than multiple rings. 

In fact, if the 13.8 day period is due to disk accretion, 
one can loosely estimate the lifetime of the accretion disk 
from the epoch- folding search results. In the epoch- folding 
searches of the first 1780 days of the BATSE data, the 
width of the x 2 peak around 13.8 days indicates that the 
peak originates in a single periodicity. 

However, in the case that a series of short-term transient 
disks formed and disappeared successively in the given 
epoch, the long term periodicity is likely made of multiple 
periods. This results because the long term periodicity 
caused by the precession of each disk could be slightly dif- 
ferent from one another, although it could be roughly the 
same as long as each disk forms and operates in similar 
conditions. Furthermore, it is not guaranteed that the pe- 
riodicity in X-rays from each disk would be in phase unless 
the disk formation procedure is periodic. The different pe- 
riodicities and their distinct phases induced by each disk 
can contribute to wash out the modulation of the folded 
light curve in epoch folding searches, and these may have 
caused the observed weakening and spreading of the long 
term periodicity in the later part of the data. 

Consequently, a single, relatively strong 13.8 day peri- 
odicity in the early BATSE data is likely due to one or 
only a few accretion disks, and the duration of the strong 
13.8 day periodicity is of order the lifetime of an accretion 
disk in the system. It is not easy to estimate the exact 
duration of the 13.8 day periodicity in the data. However, 
judging from the significance of the x 2 over ~ 300 days 
or 1700 days of data, we estimate that the lifetime of the 
transient disk could be larger than a few 100 days. 

One can also estimate the disk lifetime in relation to the 
time scale T v ; sc of viscous damping of the disk. 

T visc ~ 3.5 days a^M^R^M 1 ^, 
where a is the viscosity parameter, Mx the neutron star 
mass in M©, M\q the accretion rate in 10 16 g/sec, and 
i?io the radial distance from the neutron star in 10 10 cm 
(Frank, King, & Raine 1992). 

The accretion rate can be calculated from the observed 
flux Fx ~ 10 36 erg/sec, and Fx is mostly due to stellar 
wind capture. Assuming that the radiation efficiency e is 
~ 0.1 and that ~ 14% of Fx is from disk accretion, we get 
Mi6 ~ 0.16. The range of R for the disk can be estimated, 
too. The disc is cut off at the tidal radius i?tidos where the 
tidal torque first changes sign, and thus the outer radius of 
the disk r Q — -Rtidos- According to Wijers & Pringle (1999), 
-Rtidos ~ 0.87 Rl, where Rl is the volume-averaged Roche 
radius of the accretor. Eggleton (1983) calculated Rl in 
terms of the binary separation r& and the mass ratio qx (= 
Mx/Mo = 1/?, where Mo is the mass of the companion 
star). Assuming r b w 31i? Q and q « 12 for 4U 1700-377 
(Rubin et al. 1996), we get r a ~ 3.8 x 10 11 cm. Finally, 
according to Rubin et al. (1996), M x « 2.6 M©. 

Now we can estimate the disk lifetime from T v i sc at R 
= r Q . In 4U 1700-377, at R = r a , we gets T visc - 1900 
- 4600 days for a typical a in the range of 0.1 - 0.3. As 
mentioned earlier, the inward drift in the disk is mainly 
driven by the viscosity, and therefore, such a long viscous 
time scale implies that the disk can have a long lifetime 
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Fig. 8. — The intensity distribution of the first 1780 days of the BATSE data (the solid line) and the next 1541 days of the data (the 
dotted line): when the 13.8 day periodicity is present (the first 1780 days of the data), there are more occurrences of high count episodes. 
The average flux in the given unit is 0.0581 for the first 1780 days of the data and 0.0498 for the next 1541 days of the days. 



(> 1000 days). 

In summary, 4U 1700-377 can have a transient disk that 
can last > 1000 days due to frequent Roche-lobe over- 
filling (<~ 3.4 days). If the precession of such an accre- 
ton disk induces the 13.8 day period in X-ray emission, 
the change of the periodicity is likely related to the sta- 
bility of the precession. In addition, in this picture of 
a frequent Roche-lobe overfilled accretion disk model, in- 
stability of the accretion disk itself can contribute to the 
"weakening of the 13.8 day periodicity in X-ray emission. 
For instance, the system may have transitioned into a less 
frequent Roche-lobe overfilling state that allows the for- 
mation of very short-term transient disks only. This is 
consistent with the observed weakening and spreading of 
the long term periodicity in the later part of the data. 

3.2. Radiation driven disk warping and precession 

In the radiation-driven, warped-disk precession model, 
the radiation from the central source provides a torque to 
warp the disk and subsequently allows the precession of 
the disk. The instability criterion for the radiation-driven 
disk warping is 

167r 2 ?7 2 r b 



n 



> 



GMx/c 2 ~ e 2 To 

where n — v<ilv\ is the ratio of the two effective viscosi- 
ties associated with the shear parallel (v\) and vertical 
(V2) to the disk plane (OD). Under the assumptions of the 
model by OD, r\ can be expressed as a function of the vis- 
cosity parameter a. In Papaloizou & Pringle (1977), the 
ratio, r /rh, can be given by a function of the mass ratio 
q. Therefore, once a and e are known, one can determine 
the instability of disk warping from r& for a given q. Fig- 
ure 9, taken from OD, shows the instability condition in 
rb versus q (rb is now in GMx/c 2 ). 

In the figure, OD pointed out that systems near the in- 
stability limit tend to show a steady precession (e.g. Her 
X-l, SS 433, and LMC X-4). The systems that lie away 
from the instability limit exhibit somewhat unstable pre- 
cession (e.g. Cen X-3, SMC X-l, 4U1907+09). 

In the case of HMXRBs (q ~ 10), the instability con- 
dition becomes n,/10 6 > 3.5 for a = 0.3 and e = 0.1 
(OD). In the case of 4U 1700-377, assuming r b w 31 R & 



and M x « 2.6M Q (Rubin et al. 1996), r 6 /10 6 is - 5.4. 
Or if r b w 36 R Q and M x w 1.8 M Q (Heap & Corcoran 
1992), r fc /10 6 is - 9.4. While 4U 1700-377 does satisfy 
the instability condition of disk warping, it lies somewhat 
away from the instability limit. Therefore, according to 
OD, 4U 1700-377 belongs to the group of unstably pre- 
cessing HMXRBs such as Cen X-3, SMC X-l, and 4U 
1907+09. Therefore, if the 13.8 day period of X-ray flux 
from 4U 1700-377 is due to a radiation driven disk preces- 
sion, the instability of the periodicity is not surprising. 

It should be noted that in the above comparison with 
the model, many parameters assumed for 4U 1700-377 are 
not well constrained. In particular, a, a critical parame- 
ter of the model, is unknown for the disk of 4U 1700-377. 
Thus, as in OD, we assumed a — 0.3, which is a typi- 
cal value for a nominal disk and this assumption allows a 
direct comparison with the model described in OD. 

The precession period estimate given by OD also indi- 
rectly supports the idea of the unstably precessing disk 
in 4U 1700-377. The estimate in OD is given only for 
steadily precessing systems and the estimated period Pi on g 
depends mostly on Mx and M. For q > 1, the range of 
precession frequency w p is 

w p ~ (3.2 - 7.0) x 10- 7 C 7 2 /W (M/2Tr) 3 / 10 (GM x )- 3/2 c 5 / 2 

where Cj is a constant relating the second vertical mo- 
ment of the density and the surface density of the disk 
(Cx — 4.4 in cgs units ; refer to Eq. 9 in OD). The pre- 
dicted precession period given for Her X-l in OD allows 
one to estimate the precession period for other steadily 
precessing system. 

/ 40 \ 3/1 ° / M- ^ 3/2 
P long = (22 - 47 days)U— ' 

V-M16/ 



x 
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If 4U 1700-377 is a steadily precessing system, the esti- 
mated range of the long term period is ~ 290 - 620 days 
for M x = 2.6 Mq and M 16 ~ 0.16. Therefore, the ob- 
served 13.8 day period is too short for 4U 1700-377 to be 
a steadily precessing system according to the model. How- 
ever, OD also stated that the prediction is subject to large 
errors due to large uncertainties in the model parameters 




Fig. 9. — Stability of X-ray binaries to radiation-driven warping (a = 0.3 and e = 0.1), taken from Ogilvie & Dubus (2001). The lines 
represent the low boundary values of r;, for the instability condition for a given q, depending on the bending modes. Two estimates for 
4U 1700-377 in the text are marked with an X. HMXRBs arc shown with circles, soft X-ray transients with diamonds, other low mass X-ray 
binaries as squares. Black hole candidates appear as open symbols and neutron star primaries as solid symbols. Crosses represent systems 
for which a long term period has been reported. See Tabic 1 in Ogilvie & Dubus (2001) for the list of the binaries in this figure. 



and for 4U 1700-377, M 16 is somewhat poorly constrained 
here. 

3.3. Magnetically driven disk warping and precession 

Lai (1999) introduced a magnetically driven warp- 
ing/precession model for superorbital periods in X-ray bi- 
naries. In this model, strong magnetic fields in the magne- 
tosphere, in particular ones penetrating the accretion disk, 
exert a torque necessary for disk warping and precession. 
In magnetically driven disk warping an inner portion of the 
disk will be warped and as a result the disk will precess. 

In the case of 4U 1700-377, persistent pulsation has not 
been observed. Therefore, the existence of a strong mag- 
netic field in the system is not confirmed. However, the 
hard X-ray spectrum implies that the system may have 
strong magnetic fields. For example, a hard tail in the 
hard X-ray spectrum is the characteristic of a massive bi- 
nary system with an X-ray pulsar (Laurent et al. 1992), 
and the BeppoSAX observation showed some hints of a 
broad cyclotron absorption feature (Reynolds et al. 1999). 

The instability condition of magnetically driven warp- 
ing is that the magnetic warping gross rate T w should be 
larger than the viscous damping rate (Lai 1999), or 

r- 1 < r> 2 = T visc (r), 
where r < r m and r m is the radius of the magnetosphere 
boundary, and 

r m = (3.4 x 10 8 cm) r] n% 7 M'^' 7 M~ 2 ' 7 , 

where r\ is a constant of order unity and /X30 the magnetic 
dipole moment of the neutron star in 10 30 G cm 3 (Lai 
1999). 

For example, in the case of 4U 1700-377, T visc ~ 6 day 
at r ~ 0.01 r Q . Since T" 1 oc r 49 / 10 and T visc cx r 5 / 4 , if 
r^ 1 < 6 day at r = 0.01 r Q , T" 1 < T visc (r) at r < 0.01 r a . 
Then, the magnetic fields can warp the disk at r < 0.01 r a . 
For 4U 1700-377, assuming M XA = 1.86 (= M X /IAM Q ), 

Myj ~ 0.016, and rj ~ 1, we get r m ~ (l.OlxlO 9 cm) n% 7 . 
Therefore, in order to have magnetically driven warping 
at r < 0.01 r Ql ^ 30 > 10 from r m > 0.01 r D . 



Lai (1999) provides T w and the magnetically driven pre- 
cession frequency v prcc in terms of various parameters such 
as the magnetic dipole moment. Because of the lack of the 
confirmed persistent pulsation from 4U 1700-377, it is dif- 
ficult to estimate basic information related to the magnetic 
field of the system such as the magnetic dipole moment. 
Thus, it is even more difficult to estimate T w and f proc . 
However, we can get a magnetic-field independent rela- 
tion from the ratio of v vrcc and T w (Eq 6.3 and 6.13 in Lai 
(1999)). 

where /3 the angle between the neutron star spin axis 
and the disk rotation axis, and 9 the angle between the 
neutron star spin axis and its magnetic moment, D = 
max(^ V 2 / — 1, ^/2H/r m ) (H is the half-thickness of 
the disk and £ specifies the azimuthal pitch of the field 
line (in general \(\ < 1). 

For a precession of the whole disk, a modal analy- 
sis is required for a proper evaluation of cos/3. In gen- 
eral, /3 will alternate between its maximum and minumum 
along the azimuthal direction. Thus, azimuthally averaged 
cos/3 < 1. Assuming £ ~ 1, D ~ 0.2 (Lai 1999), cos/3 - 1 
and J^proc = 1/13.8 days, the above relation becomes 

T^ 1 - (14.1 day) tan 2 9. 

For a = 0.3, M 16 = 0.16, and M x = 2.6M , 

T visc ~ (20.2 day) r% A ~ (1.15 day) ^% 7 , 

at r — r mi where r and r m are in 10 10 cm. Therefore, if 
a magnetically field driven disk precession caused the 13.8 
day period in the X-ray emission from 4U 1700-377, 

tan 2 9 < 1.43 r^ /4 - 0.0815 . 

This result implies that the weaker the field is, the more 
likely aligned the neutron star spin axis and the mag- 
netic dipole axis are. For example, if 4U 1700-377 de- 
velops magnetically driven disk warping at r < 0.01 r a 
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(T" 1 < 6 day), the above relation reduces to 9 < 33° and 
fiso ~ 10. Or in the case that relatively weak magnetic 
fields are present, the warping would occur in a smaller 
portion of the disk. So, if the warping occurs at r < 0.002 
r a , T^ 1 < 0.81 day, 9 < 13°, and ^ 30 ~ 0.61. If there 
is indeed such an alignment, it explains the lack of the 
observed persisent pulsation from 4U 1700-377. 

In disk precession models, the central source will be suc- 
cesively seen in 'on' and 'off' states due to the varying 
obscuration by the disk. As a result, a large modula- 
tion is expected in the detected light curve, as seen in 
the well-established disk-precession sources like Her X-l, 
LMC X-4, and SMC X-l. However, 4U 1700-377 exhibits 
a somewhat small modulation in the folded light curve 
at the observed long term period (Figure 6). The small 
modulation can be explained consistently in the disk pre- 
cession models. First, there can be some extended emis- 
sion regions such as disk corona sources (Ogilvie & Dubus 
2001). Second, in 4U 1700-377, it is possible that due to 
the inherent instability of disk precession (or the disk it- 
self), for a certain fraction of precession cycles the central 
source is not completely hidden by the disk even at the 
maximum obscuration. An alternative explanation is the 
variable-accretion rate model, which is favored for the ob- 
served long term period in LMC X-3 (Wilms et al 2001). 
Further studies such as the spectral variation analysis are 



required to assess whether such a model is compatible with 
4U 1700-377. 

4. CONCLUSION 

We studied the long term - 13.8 day periodicity in X-ray 
flux from 4U 1700-377. The epoch folding analysis indi- 
cates that the periodicity was relatively strong and con- 
sisted of a single period in the early data of the BATSE 
observation. The long periodicity feature weakened and 
spread over multiple periods after ~ JD 2450000 when 
RXTE started its observation. 

We investigated the origin of the 13.8 day period in X- 
rays using disk precession models. Disk precession models 
are a natural way to explain the 13.8 day period and its 
instability. In particular, the instability of the 13.8 day 
period is consistent with the radiation-driven warped-disk 
precession model. Alternatively, the magnetically-driven 
warped-disk precession model can provide a constraint on 
the possible alignment of neutron star spin axis and the 
magnetic dipole axis, which can explain the lack of the 
observed pulsation from the system. 
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